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QUOTE : 

« Quels que soient les progrès des 
connaissances humaines, il y aura 
toujours place pour l'ignorance et 
par suite pour le hasard et la 
probabilité. » 
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Résumé : 


Lorsque des observations changent de comportement à 
partir d'un certain moment, il peut être utile de connaitre avec 
précision le moment où se produit ce changement. L'objet de 
cet article est de présenter d'un point de vue théorique 
différents modèles probabilistes de rupture puis d'illustrer 
l'utilisation de ces modèles par deux exemples. Le premier 
concerne la détermination du seuil d'élasticité d'un matériau, de 
son module de Young et des paramètres du comportement 
plastique grâce à des données expérimentales obtenues lors 
d'un test uni axial de traction. Pour certains matériaux, la 
courbe de traction ne présente pas de rupture de pente, ce qui 
rend l'identification du domaine élastique peu précis si оп 
n'utilise par un test de rupture de comportement. Le second 
exemple utilise des données bancaires et permet de mettre en 
évidence une corrélation entre un changement de la moyenne 
d'un taux d'échange et des élections. Cet exemple n'est pas 
issu de la mécanique, mais il y a certainement des 
phénomènes liés à la mécanique qui peuvent relever de ce 
type de rupture. 








І. INTRODUCTION : 


En science des matériaux, la rupture ou fracture d'un matériau est la 
séparation, partielle (comme une crique ou une fissure ou une brisure) 
ou complète, en deux ou plusieurs pièces sous l'action d'une contrainte. 


Une rupture peut être souhaitée par le concepteur de la pièce comme 
dans le cas de la conception de dispositifs de sécurité ou au contraire 
celui-ci cherche à éviter cette rupture en mettant en adéquation la 
fonction de cette pièce avec les dimensionnement et choix des 
matériaux utilisés et des procédés de réalisation. 


А côté de la méthode des moindres carrés ordinaires (MCO), la méthode 
du maximum de угаіѕетЫапсе/М\ permet aussi d'estimer les 
paramètres d'un modèle de régression, sous l'hypothèse que la vraie 
(loi) distribution desdits paramètres est connue. Si le principe pour les 
MCO est de trouver le paramètre qui minimise la somme des carrés des 
erreurs, la méthode du maximum de vraisemblance cherche par contre à 
trouver le paramètre à même (ayant une forte probabilité) de reproduire 
les vraies valeurs de l'échantillon (celles réellement observées) 


La Méthode MV va nous permettre d'estimer les paramètres du chaque 
modèle de rupture. 


Il. PRENSENTATION DE TROIS MODELES DE RUPTURE : 


Les données observées constituent une suite de variables aléatoires Х,...., Х, qui 
sont susceptibles de changer de loi après les k premières observations. 


Nous considérons le problème de test suivant 


— hypothèse nulle : les observations X,...., X, sont indépendantes et identiquement 
distribuées suivant une loi F, 


- hypothèse alternative : il existe un instant k , compris strictement entre 1 et n (ap- 
pelé instant de rupture ) tel que : les observations X4, Х,..., X; sont indépendantes 
et identiquement distribuées suivant une loi Py, et les observations X,......X, sont 
indépendantes et identiquement distribuées suivant une loi P;,. 


Les paramètres du modèle sont estimés à partir de la fonction de vraisemblance as- 
sociée à chaqu'une des hypothèses (méthode du maximum de vraisemblance). Sous 
l'hypothèse alternative, les estimateurs des paramètres 0; еї 0, vont dépendre de l'ins- 
tant k. Autrement dit, à chaque instant k nous avons une vraisemblance, donc pour 
estimer le paramètre k, nous considérons la fonction appelée rapport des vraisem- 
blances maximales (ou processus des vraisemblances). 
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Les paramètres du modèle sont estimés à partir de la fonction de vraisemblance as- 
sociée а chaqu'une des hypothèses (méthode du maximum de vraisemblance). Sous 
l'hypothèse alternative, les estimateurs des paramètres 0, et 4, vont dépendre de l'ins- 
tant k. Autrement dit, à chaque instant k nous avons une vraisemblance, donc pour 
estimer le paramètre k, nous considérons la fonction appelée rapport des vraisem- 
blances maximales (ou processus des vraisemblances). 


Dans ce travail nous supposons que les observations dont nous disposons suivent des 
lois normales. Ce choix est arbitraire, mais d'une рагі, pour la loi normale, le processus 
des vraisemblances est simple à exploiter et d'autre part, la loi normale est assez 
pertinente pour modéliser des fluctuations liées à des erreurs de mesure 

Nous désignons par ‹,, „г la densité de la loi normale de paramètres џ et о? (notée 
Л (р, о?2)). Cette densité est définie par : 








(x) = 1 f{ (2-р)? 
Фрс?) = VE exp | 202 ` 


l'hypothèse nulle sera désignée par Ho et l'hypothèse alternative sera désignée 
par Hı. 


1 Les Modèles 
Nous allons nous intéresser à trois modèles. Le premier est utile pour détecter la rup- 
ture correspondant à un changement de moyenne, le deuxième pour détecter un chan- 
gement de moyenne et de variance et le troisième est utile lorsque la moyenne reste iñ- 


changé et que seule la variance change. L'un des domaines d'application du troisième 
modèle est la modélisation des phénomènes climatologiques. 


Modele de detection de la moyenne Се modèle sera appelé imodèle M1 


Ó. Ë tont ana Ta Аи 1 
Ho: Х\,...,Х„ suivent une loi N (Ho, 05) 


(ML Advert suivent une loi N (pn, 01) 


r : эи, р А Е | _ у 
X p11...., Xn suivent une loi N(p., a?) 


Modèle de détection de la moyenne et variance Се modèle sera appelé modèle М2 


— (sj 1 


П: X1,...,Xn suivent une loi N (po, ai) 
(M2) Ху... suivent une loi №(ру, 02) 


1: Хрл... Ха suivent une loi N (2,03) 


Modèle de détection de la variance Се modèle sera appelé modèle АЗ 


Не: X1,...,X, suivent une loi N (to, ad) 
(M3) <X1...., Xk suivent une loi N (pm. ai) 


1 : Хь... Xn suivent une loi N'(p, 05) 


L'estimation des paramètres de chaque modèle est fait par la méthode du maximum 
de vraisemblance. 


2 Estimation des paramètres du modèle M1 
Sous l'hypothèse Ho la fonction de vraisemblance Lo(uo. оо. т) ой x = (z1..... та), est 


définie par : 


n 1 | 1 n | 
Lg = Y un 92 (Tj) = == >>> ехр | ==> т: — H 1 
La Це ( i) (V2r}" (сп? р | 202 2 (v Ho)” `) ) 
Afin de pouvoir de calculer le maximum de la fonction de vraisemblance par rapport à 
цо et сг, on définit la fonction "Log-Vraisemblance” par : £ = Log L. On а: 


olko, 06, z) = —nLog(V2r) – = Log(a8) - z Dle- п)? 2) 
En dérivant par rapport à ро et т; on obtient : 


| Оѓо _ 1 
dua => De Е Ий) 
QÉ n | 
=> = ad D (vi - i = pta)? 
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Le système 3) admet la solution : 


lo = — >» X; = An 
4) 


En reportant 4) dans 1) on obtient l'expresssion de la vraisemblance maxi- 


male : 
еее "2 
LX) = Цо, 63 X) = rs 5) 
Sous l'hypothèse H, la fonction de vraisemblance 2, (р, р, е. k, x) est définie par : 
[a її 
= IEA (Ti) I T po œ (Ti) 
i=1 i=k+1 
1 рн 
AUCUN (= 201 Эб mr) °) 
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1 1 re 
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1 | 
_ БЕ TERA ехр = 707 (Le. — ju) + Уе = ро} | 


et l'expression de la fonction Log-vraisemblance f(m, Hao, Е, x) : 


f, = -nLog(v2r) - < z Logo = = = (Dei - p)” + У (=; — mx) 7) 


i=] i=k+l 


d'où 


En dérivant par rapport à y, , р et с? on obtient : 


tte ий 


dun О] i=l 
Ой, 1 ы 
с === (т; — и») 8 
ðp oi an ) 
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k 
ыт (SE - m) + $ e-m) 
Le système 8) admet la solution : 


| lÊ __ 


№ i=l 


e 1 TÏ => 
lir(k) = => Уу n-i 9) 


[=K аар 


1 
a =È (ÈX; - Xa)?+ È (X;- Xp) 
En reportant 9) dans 6) on obtient l'expresssion de la vraisemblance maxi- 
male : 
,-n/2 
Lilk, X ) = -——=———- 10) 
OO (МЭл)" (i(k) 


— f;) 


D'où le rapport des vraisemblances maximales (processus de vraisemblance) pour le 
modèle ,M1 : 
Li(k X) _ sé 

АХ) 





i(k, X) = ч 





ilk) 


L'estimateur de k est déterminé par : k = Arg Max Ri(k, X). 


3 Estimation des paramètres du modèle H2 


Sous l'hypothèse Н, la fonction de vraisemblance 2, (ш, ив. 07,03, k, x) est définie par : 


k п 
= П Fu ат (Ti) П P pia т (ri) 


i=l i=k4l 


RQ ©” (-74 © mr) 12) 
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D'où l'expression de la fonction Log-vraisemblance ѓ, (ш. ш. от. 03, k, т): 


н. 
& = -nLog(V2r) – = © Log(o 1) - —— 09(02) 
13) 
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= УЫ р) — Ез ) , (zi — pa)” 


i=] i=k+H+1 
En dérivant par rapport à џ, , ш, ст et a$ on obtient : 


at | 
Эш (ш, из, 02, 02, k, £) = 
at | 

e (H, Has 1; T3, k, т) = 
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= дот, H2, Ti, ,03, k, z) = = 


Le systeme 14) admet la solution : 
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| E oS: e - 
fir(k) = T D Х; = Xx 
i= 1 
1 п r 
Ba) — == у; X-X, 
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, 15) 
| 1 ò - 
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Te > = 
En emportant 15) dans 13) on obtient l'expresssion de la vraisemblance 
maximale : 
e-n/2 
Lilk, X) = 16) 


(V2r)" (67(k))"? (63(k))">^ 
D'oü le rapport des vraisemblances maximales (processus de vraisemblance) pour le 
modèle M2 : 





кє ne ЧЕХ) l'E PE 
Ri(k, X) = 21555 = Е D РА 


L'estimateur de k est déterminé раг : К = Arg Max Ri(k. X). 
< кеп 


4 Estimation des paramètres аи modèle M3 


Contrairement aux deux modeles précédents les expressions des estimateurs de ce 
modèle ne sont pas explicites. д. (k), 22(&) et 4; (k) vérifient les relations suivantes : 








1 k 1 D: | | 
ï (k) = = X; - р k “== Xi — À k = 0 
fir (k) ETS 2 ( lir (k)) 52(F) 2% iir (k)) 
à RS sg 
ilk) = T >, (X; — lk) 18) 
i 1 = ka ; 2 
82484 = >, (Xi- fnu(k)) 
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COMMENTAIRE : 


& Tous ce calcul n'a pour but que pour nous donner une estimation des 
Paramètres du chaque Modèle, ça nous donne aussi une application du 
Méthode du maximum du vraisemblance. 

& Pour la plupart des lois de probabilité usuelles, l'estimateur du maximum de 
vraisemblance est défini de façon unique, et se calcule explicitement. Sur le 
plan théorique, il présente de nombreux avantages. Sous des hypothèses 
vérifiées par de nombreux modèles courants. 

& La méthode du maximum de vraisemblance est donc théoriquement la 
meilleure des méthodes d'estimation. 


Ili. APPLICATION : 


1. GENIE CIVIL: 


Un problème récurrent lors de l'étude du comportement plastique d'un matériau est 
la détermination précise du domaine élastique et du module de Young. Pour certains 
matériaux ne présentant pas de rupture de pente lors d'un test uniaxial de traction, 
cette détermination n'est toujours facile. Lutilisation d'un test de rupture s'avère alors 
utile. 

La variable aléatoire gaussienne X; représente alors l'écart entre la contrainte me- 
surée с; et la contrainte évaluée en fonction de la déformation (=;) et selon le modèle 
considéré : élastique linéaire ou plastique. La déformation est supposée croitre avec 


l'indice i. Pour les à avant k on cherche donc une fonction linéaire с(=;) = а=; + 8 
où о va donner le module de Young et où 3 peut correspondre à un léger décalage 
du à l'appareil de mesure. Pour les i après k on cherche une fonction non linéaire de 
gi- Si on choisit pour décrire ce comportement non linéaire un polynôme du second 
degré, on obtient un résultat peu satisfaisant : le point de rupture trouvé est représenté 
sur la figure 5.1 par . Par contre si on choisit une fonction correspondant mieux 
au comportement plastique observé, en choisissant par exemple une repésentation 
de la forme о(є;) = а,/©; + bei + с, plus conforme au comportement des données 
expérimentales, on obtient comme point de rupture le point Е sur la figure 5.1. Les 
coefficients a, 3, a, b, с, sont obtenus simultanément. Comme pour chaque valeur de k 
il faut calculer les coefficients des modèles, pour des raisons de temps de calcul il faut 
approché les modèles par une approximation qui soit linéaire par rapport aux coeffi- 
cients afin de pouvoir utiliser pour chaque valeur de К les moindres carrés linéaires. 
courbe de traction : inconel 
1200 
РЕА. 
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2. ANALYSE DU DONNEES : 


Pour cette deuxième application, même si elle est loin du 
physique, elle est intéressante car elle nous donne un autre domaine 
d'application du cette théorie probabiliste. 


Dans cette partie, nous présentons une application pratique de l'analyse aux données 
concernant le taux d'échange Franc-Deutsche Mark (qui n'existe plus) durant l'année 
1993. La série utilisée contient 256 observations (du 1 Janvier 1993 au 24 Décembre 
1993) d'une moyenne égale à ĝo = 1, 3472 et d'un écart-type égale а 4/22 = 1.4473. 
l'objet de cet étude est d'analyser l'influence des élections françaises sur le taux 
d'échange Franc-Deutsche Mark. 

La figure suivante repésente les données brutes par ordre chronologique. 

Pour gommer les petites fluctuations et mieux visualiser un point de rupture nous 





- Présentation graphique des données brutes 


avons regrouper les données en 13 et 20 classes. On constate qu'indépendamment 
du nombre de classes les deux histogrammes gardent la même allure et présentent un 
point de rupture que la méthode ci-dessus va permettre d'identifier. 
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- Répartition des données en 13 classes — Répartition des données en 20 classes 
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La figure représente le rapport des vraisemblances maximales donné par la rela- 
tion 


R, (K) 





- Fonction de vraisemblance 


On remarque que le maximum de la fonction de vraisemblance est atteint pour la 61°" 
observation, ce qui correspond à la date 26 Mars 1993. La moyenne de la série avant 
cette étant égale à 3.4254 еї celle d'après étant égale а 0, 6971. 11 est intéressant de 
constater que la date 26 Mars 1993 coïncide avec les légistatives de 1993. 


COMMENTAIRE : 


Personnellement je trouve ce modèle probabiliste très efficace et utile mais il existe 
aussi d'autre méthodes récentes et plus précis. 


Solution Proposé : 


Nous proposons de mettre en œuvre numériquement les approches statiques par 
l'intérieur et cinématique par l'extérieur du calcul à la rupture à l'aide d'éléments finis 
dédiés pour des structures de plaque en flexion et de coque en interaction 
membrane-flexion. Le problème d'optimisation correspondant est ensuite résolu à 
l'aide du développement, relativement récent, de solveurs de programmation conique 
particulièrement efficaces. 


Futur du Projet d'Optimisation : 


L'Optimisation joue un rôle très Important pour chaque individu dans tous les 
domaines. On pratique cette technique quotidiennement sans savoir, je cite 
l'exemple d'un élevé X à la semaine du contrôle, il est demandé faire un emploi pour 
optimiser sans travail pour obtenir un maximum (une très bonne note). 
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Au niveau des entreprises européennes, elles adoptent une stratégie très efficace 
schématisé comme suit : 


Créer des conditions du travail excellente Un employé satisfait Э Qualité du 
produit élevé Оп client satisfait > Entreprise gagne Plus. 


Récemment Espagne voulait essayer une stratégie qu'elle s’agit de donner aux tous 
les employés un week-end prolongé (Vendredi de plus), tous cela pour assurer le 
confort pour son peuple pour maximiser la rentabilité (win win situation). 


Apres citez ces exemples parmi les milliers qui existent, je pense et je confirme que 
les techniques d'optimisation doivent être donné plus d'importance par tout le monde 
donc si le gouvernement ou les entreprises investit dans ce domaine plus en plus ils 
vont avoir des résultats très satisfaisantes. 


IV. CONCLUSION : 


e Modèle nous permet de détecter le point de rupture donc il 
nous facilite la tâche dans plusieurs activités et diminue les 
risques et peut augmenter les gains. 


Il y a une large littérature sur la modélisation du phénomène de 
rupture en probabilité nous avons limite notre présentation à 
l'hypothèse d'une loi de probabilité normale саг cela se prête bien à 
l'estimation d'erreur de mesure mais il existe d'autres modèles bas ` 
es sur des hypothèses différentes. Il faut donc savoir дие ces 
méthodes existent et qu’elles donnent des résultats satisfaisants avec 
une procédure somme toute assez simple. 
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